have such a distribution. Disk models produce unacceptable 
quadrupole moments for distant sources or unacceptable 
(V/ Vmax) for nearby sources’*. Galactic halo distributions’* 
must be at least 50 kpc in radius without strong central condensa- 
tion to satisfy the BATSE limits on the dipole moment. Populat- 
ing a halo with neutron stars (the currently favoured candidates) 
may prove difficult. Nearby extragalactic models will have 
difficulty with the lack of correlations with M31 and the Virgo 
cluster'®. The isotropy and V/ Vmax results can be reproduced 
by cosmological models!”!*, which require isotropy but yield 
V/Vmax<0.5 because of redshift effects. These models will 
require further development to address details of the burst 
phenomenon, such as energy spectra, occurrence rate and 
luminosity. O 
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MILLISECOND radio pulsars, which are old (~10° yr), rapidly 
rotating neutron stars believed to be spun up by accretion of matter 
from their stellar companions, are usually found in binary systems 
with other degenerate stars". Using the 305-m Arecibo radiotele- 
scope to make precise timing measurements of pulses from the 
recently discovered 6.2-ms pulsar PSR1257 + 12 (ref. 2), we demon- 
strate that, rather than being associated with a stellar object, the 
pulsar is orbited by two or more planet-sized bodies. The planets 
detected so far have masses of at least 2.8 Mo and 3.4 M,, where 
Mg is the mass of the Earth. Their respective distances from the 
pulsar are 0.47 AU and 0.36 AU, and they move in almost circular 
orbits with periods of 98.2 and 66.6 days. Observations indicate 
that at least one more planet may be present in this system. The 
detection of a planetary system around a nearby (~ 500 pc), old 
neutron star, together with the recent report on a planetary com- 
panion to the pulsar PSR1829— 10 (ref. 3) raises the tantalizing 
possibility that a non-negligible fraction of neutron stars observ- 
able as radio pulsars may be orbited by planet-like bodies. 

The 6.2-ms pulsar PSR1257+ 12 (Fig. 1) was discovered dur- 
ing the search at high galactic latitudes for millisecond pulsars 
conducted in February 1990 with the 305-m Arecibo radiotele- 
scope at a frequency of 430 MHz (ref. 2). The characteristics of 
this survey and the details of data analysis are described else- 
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where*. The confirming observations made on 5 July 1990 have 
been followed by routine pulse timing measurements of the new 
pulsar. A total of 4,040 pulse time-of-arrival (TOA) observations 
have been accumulated so far, with the Arecibo radiotelescope, 
the 40-MHz, three-level correlation spectrometer and the Prin- 
ceton Mark IH pulsar processor at 430 MHz and 1,400 MHz. A 
typical uncertainty in the TOAs derived from 1-min pulse 
integrations is ~15 ps. 

The standard analysis of the timing data has been carried out 
using the model fitting program TEMPO® and the Center for 
Astrophysics Solar System ephemeris PEP740R. With a growing 
time span of the TOA measurements, it had gradually become 
clear that the TOAs showed an unusual variability superimposed 
on an annual sinusoidal pattern caused by a small (~1’) error 
in the assumed pulsar position. To separate these effects unam- 
biguously, a timing-independent, interferometric position of 
PSR1257+ 12 was measured with the Very Large Array (VLA) 
in its A-array configuration on 19 July and again on 18 September 
1991. The ~0.1” accuracy of the resulting pulsar position was 
achieved by referencing the fringe phase to a point-source cali- 
brator 1.7° away. 

A least-squares fit of a simple model, which involved the 
pulsar’s rotational period, P, and its derivative, P, as free para- 
meters and the fixed VLA position (Table 1), to the timing data 
spanning the period of 486 days resulted in post-fit residuals 
shown in Fig. 2a. The residuals, which measure the difference 
between the predicted and the actual TOAs, show a 
quasiperiodic ‘wandering’ over the entire pulsar period. A closer 
examination of this effect has revealed that it was caused by 
two strict periodicities of 66.6 and 98.2 days in the pulse arrival 
times. This is further demonstrated in Fig. 2b and c, which shows 
post-fit residuals after fitting each of the two periods separately 
to the above data, assuming simple keplerian binary models 
involving a low-mass binary companion to the pulsar. Evidently, 
fitting for one of the assumed binary periods leaves the other 
one as a post-fit residual, implying that the pulse arrival times 
of PSR1257+12 are indeed affected by two independent peri- 
odicities. Further detailed analysis has shown that the 
periodicities are independent of radio frequency and that other 
millisecond pulsars routinely observed at Arecibo with the same 
data acquisition equipment show no such effect in their timing 
residuals. 

Millisecond pulsars are extremely stable rotators. Systematic 
timing observations of objects like the 1.5-ms pulsar 1937+ 21 
(ref. 6) have not revealed any timing noise, quasiperiodic TOA 
variations or ‘glitches’ at the level often found in the population 
of younger pulsars and believed to be related to neutron star 
seismology’. The frequency independence of the amplitude of 
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FIG. 1. The average pulse profile of PSR1257 +12 at 430 MHz. The effective 
time resolution is ~12 ws. 
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TABLE 1 Parameters of the PSR1257 +12 system 


Pulsar parameters 
Rotational period, P 
Period derivative, P 
Right ascension (B1950.0, VLA) 
Right ascension (B1950.0, timing) 
Declination (B1950.0, VLA) 
Declination (B1950.0, timing) 


0.00621853193177 + 0.00000000000001 s 
1.21x1079+2x10°74 s 67+ 

12 h 57 min 33.131 s+0.015 

12 h 57 min 33.126 s+0,003 

12° 57'05.9"+ 0.4 

12° 57'06.60" +0.02 


Epoch JD 2448088.9 
Dispersion measure 10.18+0.01 pe cm? 


Flux density (430 MHz) 20+5 mJy 
Flux density (4,400 MHz) 1.0+0.2 mJy 
Surface magnetic field, B 8.8x10°G 
Characteristic age, 7, 0.8 x10° yr 


Keplerian orbital parameters 
Projected semimajor axis, & sin 1.31 +0.01 light ms 1.4140.01 light ms 
Eccentricity, e 0.022 + 0.007 0.020 + 0.006 
Epoch of periastron, To JD 2448105.3+1.0 JD 2447998.6 + 1.0 
Orbital period, P, 5751011.0 + 800.0 s 8487388.0+ 1800.0 s 


Longitude of periastron, w 252° + 20° 107°+20° 
Parameters of the planetary system 

Planet mass, M23 (Mo) 3.4/sin i 2.8/sin i 

Distance from the pulsar, d (au) 0.36 0.47 

Orbital period, P, (days) 66.6 98.2 


TOA variations in PSR1257+12 rules out propagation 
phenomena, such as those detected in eclipsing binary pulsars*”’, 
as a possible source of the observed periodicities. The integrated 
pulse profiles of PSR1257+ 12 do not show any morphological 
changes that might indicate the presence of a free precession 
of the pulsar spin axis'™!! or any magnetospheric phenomena 
at the level that could lead to periodic TOA variations of the 
measured magnitude. Consequently, the most plausible remain- 
ing alternative is that PSR1257+12 has two low-mass com- 
panions and that their orbital motion is responsible for the 
observed TOA variations. 

To analyse this exciting possibility further, we have modified 
the code of the TEMPO program to accommodate a timing 
model including multiple, noninteracting keplerian orbits. The 
result of a fit of the model including the rotational and positional 
parameters of the pulsar as well as the two, five-parameter 
keplerian orbits to the entire data set is shown in Fig. 2d. The 
post-fit r.m.s. residual of the resultant timing model is ~18 ws 
(comparable to the individual TOA uncertainties) and the 
remaining residuals contain very little systematic variation. The 
model parameters of the pulsar and its assumed two companions 
are listed in Table 1. The derived parameters have been obtained 
from standard considerations involving very low-mass objects 
in keplerian orbits around a 1.4Mọo neutron star. The pulsar 
period variations predicted by this model are displayed in Fig. 
3 together with the observed period values derived from timing 
data. Note that the presence of the planet-sized bodies orbiting 
the pulsar results in apparent period variations of only +15 ps. 
This is caused by the pulsar’s spatial motions which are charac- 
terized by the projected maximum velocity and displacement 
amplitudes of +0.7 ms”! and +900 km, respectively. 

The high quality of the two-companion model fit to the pulsar 
timing data covering several orbital cycles provides compelling 
evidence that PSR1257 + 12 possesses a planetary system consist- 
ing of at least two planet-sized bodies. The possibility of more 
planets around PSR1257+ 12 is indicated by a 0.7” discrepancy 
between the VLA and the timing positions, which is considerably 
greater than the conservative error estimates (Table 1). This 
discrepancy may arise from a bias in the timing position caused 
by a presence of a third planet with orbital period close to one 
year. Because the currently measured value of the second period 
derivative, P, of the pulsar is not significant, the effect of any 
outer planets that could be present in the 1257+ 12 system would 
be entirely absorbed in P. 

The characteristics of the 1257+12 planets are not unlike 
those of the inner Solar System. Both planets circle the pulsar 
at distances similar to that of Mercury in its orbit around the 
Sun. Assuming a random distribution of the inclinations, i, of 
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orbital planes, there is a 50% probability that i= 60°, so that 
the median values of the permissible masses are 3.2M, and 
3.9M,, respectively (see Table 1). Interestingly, the ratio of 
orbital periods, 1.476+0.001, is close to a 3/2 orbital resonance 
of the type often encountered in the Solar System, either between 
planetary satellites, between Jupiter and some asteroids, or even 
Neptune and Pluto’. Also, the similarity of the measured eccen- 
tricities combined with the ~ 180° separation of the pericentres 
of the orbits can be easily understood in terms of secular per- 
turbations of the orbital elements of the two planets (see, for 
example ref. 13). Finally, a hypothetical energy flux at the 
planets’ distance from the pulsar (~0.4 AU) due to the neutron 
star’s rotational energy loss (E = Iwo, where I is the moment 
of inertia and w =27/P) is ~4x 10’ erg cm~? s™*. This is ~30 
times the solar constant and corresponds to a black-body tem- 
perature of ~670 K, which is similar to the measured dayside 
surface temperature of Mercury (see, for example, ref. 14). 
The possibility of the existence of planet-sized bodies orbiting 
neutron stars has been contemplated in the past'*'°. The most 
recent evidence has been presented by Bailes et aP, who have 
detected a ~6-month periodicity in the timing residuals of a 
relatively young pulsar, PSR1829—10, which is interpreted in 
terms of the orbital motion of a ~10 Mz, companion. A number 
of mechanisms have been proposed to explain a planet around 
PSR1829 — 10 (refs 17-21). In the case of PSR1257+ 12, its old 
age and a low surface magnetic field (Table 1) are typical of 
neutron stars which are believed to evolve in low-mass binary 
systems and are spun up to the observed millisecond periods 
by accretion of matter from their stellar companions’. As it is 
not very likely that any primordial planets would survive this 
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FIG. 2 The post-fit residuals of pulse arrival times from PSR1257 +12. a, 
Fit for rotational parameters only (VLA position of the pulsar fixed); b, fit 
for a 98.2-day keplerian orbit (leaves a 66.6-day periodicity as residual); c, 
fit for a 66.6-day keplerian orbit (leaves a 98.2-day periodicity as residual); 
d, fit for all parameters of a-c. 
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FIG. 3 Period variations of PSR1257 +12. Each period measurement is 
based on observations made on at least two consecutive days. The solid 
line denotes changes in period predicted by a two-planet model of the 
1257 +12 system. 


kind of evolution”’, the observed 1257+12 system probably 
consists of ‘second generation’ planets created at or after the 
end of the pulsar’s binary history. Because a second supernova 
explosion is not expected to occur in a low-mass binary, such 
planets could form a stable system of bodies orbiting an old 
neutron star. Another important evolutionary constraint is pro- 
vided by the observed very low eccentricities (e ~ 0.02) of the 
planetary orbits and the near-resonance ratio of the orbital 
periods. These characteristics suggest that the planets were 
created from some form of accretion disk that would naturally 
provide means to circularize the orbits and to bring them close 
to a 3/2 resonance’”. Consequently, it seems that any plausible 
mechanism for the creation of a planetary system around a 
millisecond pulsar must provide a way to remove its stellar 
companion in a manner different from the disruption of a binary 
system caused by a supernova explosion and to retain enough 
circumpulsar matter to form planet-sized objects. 

Within the constraints provided by the observational evidence, 
it is tempting to postulate that the 1257+12 system simply 
represents one of the possible outcomes of a neutron star evol- 
ution in a low-mass binary system'. As two examples of com- 
panion-star evaporation by millisecond pulsars have already 
been detected”*, this mechanism seems to account naturally for 
the absence of such a companion to PSR1257+ 12. In this case, 
one way to supply the material for planet creation could be a 
disk formed out of the ablated stellar matter ejected from the 
binary system****. The existence of such an ‘outer disk’ around 
the eclipsing binary pulsar PSR1957+20 could explain the 
observed orbital decay of this system”>. 

The results described here strongly suggest that one of the 
nearby galactic millisecond pulsars, PSR1257+ 12, is accom- 
panied by a system of two or more planet-sized bodies. The 
pulse timing method (similar to the analysis of single-line 
spectroscopic binaries) was used to detect a +0.7ms ‘ pulsar 
‘wobble’ caused by orbital motions of the planets. In fact, 
planetary masses smaller than that of the Moon could easily be 
detected with this technique. At present, this kind of accuracy 
is entirely inaccessible to the optical methods of planetary system 
detection”. The existence of both the 1257+ 12 planetary system 
and the object orbiting the pulsar PSR1829 — 10 seems to suggest 
that planets can form under a variety of conditions. This notion 
and the possibility of a non-negligible frequency of occurrence 
of planets around neutron stars, if confirmed, will have far- 
reaching consequences for our understanding of the formation 
and evolution of planetary systems and for future strategies of 
searches for planets outside the Solar System. O 
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FOLLOWING an initial radio outburst’ which decayed on a time- 
scale of a few weeks, supernova 1987A? was undetectable at radio 
wavelengths until recently. In mid-1990, 1,200 days after the 
explosion, radio emission was once again detected? with the 
Molonglo Observatory Synthesis Telescope and the Australia 
Telescope Compact Array. Our observations since then show that 
the source has increased in strength at all monitored frequencies 
between 843 MHz and 8.6 GHz, with considerable variations in 
spectral index. The extended radio emission is centred within 
0.5 arcsec of the optical supernova, and probably lies within the 
IO n) ring imaged by the Hubble Space Telescope. We interpret 
the emission as optically thin synchrotron emission from shock- 
accelerated electrons. This is the first time the birth of a nearby 
radio supernova remnant has been witnessed, and future observa- 
tions will allow the structure of the remnant to be compared with 
the many other known radio remnants. 

Since the demise of the initial radio emission, we have 
monitored SN1987A using the Molonglo Observatory Synthesis 
Telescope* (MOST) at 843 MHz and, since mid-1989, the Aus- 
tralia Telescope Compact Array? (ATCA) at 1.5, 2.4, 4.8 and 
8.6 GHz. Renewed radio emission at the supernova position was 
first detected on 6 July 1990, using the MOST”. The flux density 
at this time was ~3 mJy. At higher frequencies, the onset of the 
radio emission was delayed, the first detection with the ATCA 
at 4.8 GHz being on 16 August 1990 (0.7 mJy). Since then, the 
MOST and the ATCA have made 70 and 40 observations, 
respectively. Each of these observations lasted ~12 hours. The 
MOST observations are affected by the strong source 30 
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